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Abstract. The endpoint of black hole evaporation is a very intriguing problem of modern 
physics. Based on Einstein-dilaton-Gauss-Bonnet four dimensional string gravity model we show 
that black holes do not disappear and should become relics at the end of the evaporation process. 
The possibility of experimental detection of such remnant black holes is investigated. If they 
0^ \ really exist, these objects could be a considerable part of the non baryonic dark matter in our 

' Universe. 
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Introduction 



Theoretical physics faces nowadays a great challenge. There is four dimensional 
^ ■ Standard Model on one side (and the additional dimensions are not required to explain 
experimental data) together with inflationary cosmology based on additional scalar 
fields Q. On the other side there is the completely supersymmetrical string/M-theory. 
Building links between those approaches is a very motivating goal of modern physics 
which could be achieved by the study of microscopic black holes. 

As General Relativity is not renormalizable, its direct standard quantization is 
impossible. To build a semiclassical gravitational theory, the usual Lagrangian should 
be generalized, which is possible in different ways. One of them is to study the action 
expansion in scalar curvature, i.e. higher order curvature corrections. At the level of 
second order, according to the perturbational approach of string theory, the most natural 
choice is the 4D curvature invariant Gauss-Bonnet term Sgb = RijkiR^ kl — ^RijR 1 ^ + R 2 

With 4D action, it is not possible to consider only Sgb because, being full derivative, 
it does not contribute to the field equations. It must be connected it with a scalar field 
to make its contribution dynamical. The following 4D effective action with second 
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order curvature corrections can be built: 



S 



d A X\J—g 



-R + 2d^(j) + X^<t>)S GB + ... 



where A is the string coupling constant. As in cosmology, the most simple generalization 
of the theory (a single additional scalar field) is not possible because while dealing with 
spherically symmetric solutions, the "no-hair" theorem restriction must be taken into 
account. 

Treating as a dilatonic field, the coupling function £(</>) is fixed from the first 
string principles and should be written exp(— 2<p) |4|, ||, which leads to : 



S 



d A xyf—g 



-R + 2d, 



+ Xe-^ScB + ... 



(1) 



Such type of actions can be considered as one of the possible middle steps between 
General Relativity and Quantum Gravity. In this paper, we show that this effective 
string gravity model and its solutions can be applied for a description of the last stages 
of primordial black holes (PBH) evaporation || |9| and suggests possible dark matter 
candidates |TIJ. This should be understood in the general framework of Gauss-Bonnet 
black hole (BH) theory g 0]. 

The paper is organized as follows: In Section II we briefly recall previously obtained 
results and point out some new features important for this study, Section III is devoted 
to the establishment of the new Hawking evaporation law (especially for the detailed 
description of last stages of Gauss-Bonnet BH evaporation), in Section IV we show 
that the direct experimental registration of such PBHs is impossible, Section V is 
devoted to PBH relics as dark matter candidates and Section VI contains discussions 
and conclusions. 



Black hole minimal mass 



2.1. Black hole minimal mass in pure EDGB model 

For the sake of completeness, main results from Ref. j7j are briefly repeated. 

Starting from the action (0), a static, spherically symmetric, asymptotically flat 
black hole solution is considered. One of the most convenient choice of metric in this 
model is 



ds 2 



Adt 2 - 



2 ° dr 2 



A 



r 2 {dd 2 



sin 



(2) 



where A = A(r), o = cr(r). 

Asymptotic expansion of the solution has the usual quasi-Schwarzschild form, 
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oo 



r \r 



ID 2 

a(r^oo) = l--— + O yr2 

Mr -> oo) = — + of- 
r \r 

where M and D are ADM (Arnowit-Dieser-Misner) mass and dilatonic charge 
respectively. Using a dedicated code, a BH type solution was obtained. This solution 
provides a regular horizon of quasi-Schwarzschild type and the asymptotic behaviour 
near this horizon is: 

A = di(r - r h ) + d 2 {r - r h ) 2 + . . . , 
a = s Q + s 1 (r-r h ) + 

= 0oo + 0i - r h ) + 2 (r - r h f + . . . , (3) 

where (r — r^) -C 1, s , 0o = e~ 2 * 00 and r/j are free independent parameters. 

After solving the equations to the first perturbation order, the following limit on 
the minimal BH size can be obtained: 



r 



inf 



h ~ VA VW60h(&a), ( 4 ) 

where A is a combination of the string coupling constants (fundamental value) and 
0/i(0oo) is dilatonic value at r^, depending upon dilatonic value at infinity which cannot 
be determined only in the framework of this model. According to this formula and 
taking into account the numerical values, the minimal BH mass has the order of Planck 
one (more precisely ~ 1.8 Planck masses 0). 

It is necessary to point out that the stability of the solution under time perturbations 
at the event horizon was described in and was studied at the singularity r s in |12 



Contribution of higher order curvature corrections was studied in 0, [16[] to show 
that, in the bosonic case with heterotic string models (the question is still open in SUSY 
II), all the new topological configurations are located inside the determinant singularity 
and, therefore, do not produce any new physical consequences. Our conclusions remain 
valid when the next higher order curvature corrections are made of pure products of 
Riemannian tensors. This topic is under additional investigation now. 

Gathering these results, it can be concluded that the solution is stable in all the 
particular points, and, therefore, at all the values of initial data set. 

2.2. Effects of moduli fields 

Generalizing the model by taking into account the effective contribution of additional 
compact dimensions in the most simple form — scalar field — the action must be given 

as 



S = J d^x^g -R + 2d ll (j)d ll (j) + 2d^d^ + ^e" 20 + 



GB 
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+ higher order curvature corrections 



(5) 



This model was studied into the details in ||13|| . For current investigations it should 
be emphasized that when the contribution of moduli field value is considered, a naked 
singularity can appear if the size of additional dimensions is greater than the BH size. 
The minimal BH mass must therefore be increased to 10 Planck masses (to avoid being 
in naked singularity region) . It is a key feature because it allows to move away from the 
Planck region and to use semiclassical approach. If these additional dimensions were 



non-compact |T4| the BH minimal mass would be much greater. 



3. Black hole evaporation law 



3.1. Probability of transition to the last stage 



According to the analysis given in Ref. |18|], the transition from prelast to last stage of 
BH evaporation is forbidden and evaporating PBHs will never reach the minimal mass 
state. The shape of the BH mass loss rate law changes and becomes the one presented 



in Fig.2, analogously to the simplified "toy model" presented in [H3|. Different types of 



similar models for BH evaporation were studied in Lovelock gravity |T^] , string inspired 
curvature expansions |30j and in many other theories. The numerical values of Gauss- 
Bonnet BH (important for experimental search analysis) will be presented in section 
D. 



3.2. Approximation to metric functions 

In the WKB approximation of the Hawking evaporation process, everything happens in 
the neighbourhood of the event horizon. As our metric functions A and a depend upon 
radial coordinate r and black hole mass M, i.e. A = A(M, r) and o = <x(M, r) (other 
variables are not important), we can use expansions @, taking into account only the 
first terms (partially neglecting the dependence upon radial coordinate r). Using [| the 
metric can be written as: 

A(M,r) = l-^ e (M) = ^-(r-2Me(M)), 

a(M,r) = ao(M). (6) 

Using the numerically calculated data, fits were performed for e(M) and <r (M). As 
we are interested mostly in the last stages of PBHs evaporation, where the difference 
from the standard Bekenstein-Hawking picture is considerable, Taylor expansions 
around M min can be used. This also helps in obtaining good fits of the metric functions 
(see Fig. 3 and Fig.4) which can be considered as polynom expansions (of M or 1/M) 
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that are valid between M = M min = 10 Mpi and M = 1000 Mpi with good accuracy. 

f MM 2 M 3 M 4 ' 1 

do = a 2 (M-M mm ) 2 -a 3 (M-M min ) 3 + a4(M-M min ) 4 -a 5 (M-M miri ) 5 , 

where (for M min = 10 Mpi) the corresponding coefficients are t\ = 10.004, e 2 = 
13.924, e 3 = 2856.3, e 4 = 25375.0, a 2 = 0.11933 * lO" 04 , a 3 = 0.30873 * lO" 07 , 
CJ4 = 0.30871 * 10- 10 , cr 5 = 0.11051 * 10~ 13 . 

Using this technique, the PBH evaporation spectra and mass loss rate were derived 
in an analytical form (valid only near the M m i n point). 



3.3. Black hole evaporation spectra in EDGB model 



In some approaches, black holes are treated as immersed in a thermal bath and the 
evaporation can be described as a WKB approximation of semiclassical tunnelling in 



a dynamical geometry. In our investigation, we follow the techniques described in [pi 
and p2| . The same method was also applied in [23fl. Some other descriptions of BH 



evaporation can be found in [24, 25 



The key idea of the method from Ref. [21] and [22] is that the energy of a particle 
changes its sign when crossing the BH horizon. So, a pair created just inside or just 
outside the horizon can become real with zero total energy after one member of the pair 
has tunnelled to the opposite side. The energy conservation plays a fundamental role: 
transitions between states with the same total energy are the only possible ones. Using 
quantum mechanical rules, it is possible to write the imaginary part of the action for 
an outgoing positive energy particle which crosses the horizon outwards from r in to r out 
as: 



Im(S) = Im 



M- 



r ut dr 
—an, 
r 



where io is the energy of the particle, H is total Hamiltonian (and total energy) and 
the metric is written so as to avoid the horizon coordinate singularity. Following 
22| , Painleve's coordinates are used. The transformation to this metric from the 



Schwarzschild one can be obtained by changing the time variable: 



t = t old + r\ 



• a* 



1 

A' 



Substituting t u into (||) one obtains 

ds 2 = -Adt 2 + 2Va 2 - Adrdt + dr 2 + r 2 dQ 2 . 



(8) 



In WKB approximation, the imaginary part of the semiclassical action Im(S), 
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describing the probability of tunnelling through the horizon is 

r out '''out Pr 

Im(S) = Im J p r dr = Im J J p' r dr, (9) 

Tin fin 

where p r is canonical momentum. 

For Gauss-Bonnet BH the radial geodesies are described by the equation 

= a, = = T<r~ V(J 2 - A- (10) 

After substituting expression ( PH ) to the equation (^) one obtains: 

f f dr r f f drduj , , 

Im 5 = Im —dH=-Im / . 11 

J J f J J a - va 2 - A 

M r in 2Me 

Substituting the expression @ extended in (0) to the equation (|TT|) , the imaginary 
part of the action can be written as: 

. 2(Af-c/)e 



/m(S) = -Im J du'i J 



o.W, ^ 2(M-u J >)e +1 < 



Changing variables with 



2(M-u)e 
Im(S) takes the form 

4(M - uj')eydy 



Im(S) = -Im I duj'i J 



y-a 



-Im J du'U(M - u)')ea J 

n V / 



dy 



r 2 1* 

M-ui' 



dJ (4(M - u/)e<rir) . 



As a result the imaginary part of the action is 
840vr 
M 2 (M -uo) 

where a is a huge expression that cannot be written here. It can be found at 
http: / / isnwww.in2p3.fr/ams/lmS.ps. 



2Im(S) = ^777 — rA2 a > 



Using the numerical values for a realistic order of M m i n around 10 Planck masses, 
the corresponding and <Tj, it is possible to find from (g) the approximate expression 
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for Im(S). As we are interested mostly in the last stages of BH evaporation where the 
influence of higher order curvature corrections is important, the limit M — M min <C 1 can 
be taken in the computations, leading to a very different spectrum than the standard 
Bekenstein-Hawking picture (where —dM/dt oc 1/M 2 ). Taking into account energy 
conservation, u can be bounded : < u < M — M min . The approximate expression of 
Im(S)(M,uj) for a given M min can then be used in the form 

Im(S) = k * (M - M min )\ (12) 

where constant k = 5 • 10 -4 in Planck unit values with a satisfying accuracy (the plot 
of Im(S) and its approximation are shown on Fig.5). 



3.4- Energy conservation and mass lost rate 



Following Ref. [22], the emission spectrum per degree of freedom can simply be written 
as: 

d 2 N Q((M — M min )c 2 — E) 

dEdt 2nh ' e Im ^ - (-1) 2 * ' 1 j 

T S (M, E) being the absorption probability for a particle of spin s and the Heavyside 
function being implemented to take into account energy conservation with a minimal 
mass M min . In this section and in the following ones, standard units are used instead 
of Planck ones as numerical results should be obtained for experimental fluxes. At this 
point, two questions have to be addressed: which kind of fields are emitted (and which 
correlative T s have to be used) and which mass range is physically interesting. To answer 
those questions, the mass loss rate is needed: 

dM r(M~M min y d 2 N E dE 
~~dt ~ Jo dEdt ' c 5 ( ' 

where the integration is carried out up to (M—M m i n )c 2 so as to ensure that the transition 
below M min is forbidden. The absorption probabilities can clearly be taken is the limit 
GME/hc 3 < 1 as we are considering the endpoint emission when the cutoff imposed 
by M m i n prevents the black hole from emitting particle with energies of the order of kT. 
Using analytical formulae ]27| and expanding exp(Jm(S')) to the first order with the 
approximation according to (|T2"D, it is easy to show that the emission of spin-1 particles, 
given by (per degree of freedom) 

dM 16 G A M Pl 



-M~<r.JFsk M(M - M ^ (15) 

dominates over s=l/2 and s=2 emission whatever the considered energy in the previously 
quoted limit. It is interesting to point out that the fermion emission around M min is not 
strongly modified by the EDGB model as, in the lowest order, exp(/m(S')) — (— l) 2s ~ 2. 
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Furthermore, if energy conservation was implemented as a simple cutoff in the Hawking 
spectrum, the opposite result would be obtained: s=f/2 particles would dominate the 
mass loss rate as the power of (ME) in the absorption probability is the smallest one. 
If we restrict ourselves to massless particles, i.e. the only ones emitted when M is 
very close to M m i n , the metric modification changes the endpoint emission nature from 
neutrinos to photons. The real mass loss rate is just twice the one given here to account 
for the electromagnetic helicity states. 

With this expression —dM/dt = f(M), it is possible to compute the mass M at any 
given time t after formation at mass Mi n n as: 

rM init dM _ 9nkh 5 c 2 1 
t ~J M 7(M) ~ 32G 4 MJ, X M^ n (M - M mm f ( } 

where only the dominant term in the limit t — > oo is taken from the analytical primitive 
of the function. As expected, the result does not depend on M init which is due to the 
fact that the time needed to go from M init to a few times M min is much less than the 
time taken to go from a few times M min to M as long as M init « 10 15 g for t ps 10 17 s. 
At time t after formation, the mass is given by: 



M w M min + 



\ SMi m G±MU 

This mass can be implemented in the emission spectrum formula: 



d*N ^32 / 8 \§ f 15 f 

dEdt^^K^) G n c 1 min 



3 



fc-3 x t^E^e 



leading to a frequency f given by 

f (M-M min y d 2 N 36 l 

If we want to investigate the possible relic emission produced now from PBHs formed 
in the early universe with small masses, this leads to a frequency around 6 ■ 10~ 18 Hz 
with a typical energy of the order of 1.8 • 1CT 6 eV. This emission rate is very small as 
it corresponds to the evaporation into photons with wavelength much bigger that the 
radius of the black hole. It should, nevertheless, be emphasized that the spectrum is 
a monophonically increasing function of energy, up to the cutoff, with a E A behaviour. 
Furthermore, it shows that, although very small in intensity, the evaporation never stops 
and leads to a mass evolution in 1/ \fi. 
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4. Experimental detection 

We investigate in this section the possibility to measure the previously given relic 
emission. Let R be the distance from the observer, z the redshift corresponding to 
the distance R, 9 the opening angle of the detector (chosen so that the corresponding 
solid angle is f2 = 1 sr), d 2 N/dEdt(E ) t) the individual differential spectrum of a black 
hole relic (BHR) at time t, p(R) the numerical BHR density taking into account the 
cosmic scale factor variations, R max the horizon in the considered energy range, t un i V 
the age of the Universe and H the Hubble parameter. The "experimental" spectrum F 
(J -1 -s -1 -sr -1 ) can be written as: 



Rmax tfN ( R 



dEdt \ c 



x , ( fl).*gw Wdfi (20) 



which leads to: 

F = -tan 2 (6) 

xk^E 4 f 
Jo 

x0 




9knh 5 c e 



(21) 

This integral can be analytically computed and takes into account both the facts 
that BHRs far away from Earth must be taken at an earlier stage of their evolution 
and that energies must be redshifted. Even assuming the highest possible density of 
BHRs (Qbhr — Qcdm ~ 0.3) and R max around the Universe radius, the resulting flux 
is extremely small: F m 1.1 • 10 7 J -1 s _1 m~ 2 sr -1 around 10~ 6 eV, nearly 20 orders 
of magnitude below the background. This closes the question about possible direct 
detection of BHRs emission. 

Another way to investigate differences between EDGB black holes particle emission 
and a pure Hawking spectrum is to study the mass region where dM/dt is maximal. 
Taking into account that the mass loss rate becomes much higher in the EDGB case 
than in the usual Hawking picture, it could have been expected that the extremely 
high energy flux was strongly enhanced. In particular, it could revive the interest in 
PBHs as candidates to solve the enigma of measured cosmic rays above the GZK cutoff. 
Nevertheless, the spectrum modification becomes important only when the mass is quite 
near to M min . Depending on the real numerical value of M min , it can vary substantially 
(increasing with increasing M min ) but remains a few Planck masses above M min . This 
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is far too small to account for a sizeable increase of the flux. The number of particles 
emitted above 10 20 eV in a pure Hawking model is of the order of 10 15 ||. This value 
should be taken with it relies on the use of leading log QCD computations 

of fragmentation functions far beyond the energies reached by colliders but the order 
of magnitude is correct. On the other hand, even if all the energy available when 
EDGB modifications becomes important were released in 10 20 eV particles (which is 
not realistic) it would generate only a few times 10 9 particles and modify by less than 
0.01% the pure Hawking flux. It would not allow to generate, as expected, a spectrum 
harder than E~ 3 . 



5. Primordial black holes as dark matter candidates 

The idea of PBH relics as a serious candidate for cold dark matter was first mentioned in 



Til • It was shown that in a Friedman universe without inflation, Planck-mass remnants 
of evaporating primordial black holes could be expected to have close to the critical 
density. Nevertheless, the study was based on the undemonstrated assumption that 
either a stable object forms with a mass around Mpi or a naked space-time singularity 
is left. Our study provides new arguments favouring massive relic objects, probably 
one order of magnitude above Planck mass and could revive the interest in such non- 
baryonic dark matter candidates. An important problem is still to be addressed in 
standard inflationary cosmology: the rather large size of the horizon at the end of 
inflation. The standard formation mechanism of PBHs requires the mass of the black 
holes to be of the order of the horizon mass at the formation time and only those created 
after inflation should be taken into account as the huge increase of the scale factor would 
extremely dilute all the ones possibly formed before. It is easy to show that under such 
assumptions the density of Planck relics is very small: 

Mpi = n PBH ?—^M}r a M:- 2 M min (22) 
a — 1 

where VLpbh is the density of PBHs not yet evaporated, a is the spectral index of the 
initial mass spectrum (=5/2 in the standard model for a radiation dominated Universe), 
M* is the initial mass of a PBH which evaporating time is the age of the Universe 
(^ 5 • 10 14 g) and Mh is the horizon mass at the end of inflation. This latter can be 
expressed as 

ilM Pl ilMp, 0.24 
8 tn 8 tpi (T RH /lMev) 2 

where Tj is the formation time and Tpu is the reheating temperature. Even with the 
highest possible value for Trh, around 10 12 GeV (according to Ref. if the reheating 



temperature is more than 10 9 GeV, BHR remnants should be present nowadays) and 



11 



the upper limit on Qpbh coming from gamma-rays, around 6- 10~ 9 the resulting density 
is extremely small: Qpi ~ 10~ 16 . 

There are, nevertheless, at least two different ways to revive the interest in PBH dark 
matter. The first one is related to relics that would be produced from an initial mass 
spectrum decreasing fast enough, so as to overcome the gamma-rays limit. The second 
one would be to have a large amount of big PBHs, between 10 15 g and 10 25 g, where 
experiments are completely blind: such black holes are too heavy to undergo Hawking 
evaporation and too light to be seen by microlensing experiments (mostly because of the 
finite size effect |[29[ ). The most natural way to produce spectra with such features is 
inflationary models with a scale, either corresponding to a change in the spectral index 
of the fluctuations power spectrum [Bof either corresponding to a step |3T] . 



6. Discussion and conclusions 

In this paper, the BH type solution of 4D effective string gravity action with higher 
order curvature corrections was applied to the description of BHRs. A corrected version 
of the evaporation law near the minimal BH mass was established. It was shown 
that the standard Bekenstein-Hawking evaporation forluma must be modified in the 
neighbourhood of the last stages. Our main conclusion is to show that contrarily to 
what is usually thought the evaporation does not end up by the emission of a few 
quanta with energy around Planck values but goes asymptotically to zero with an infinite 
characteristic time scale. 

The direct experimental registration of the products of evaporation of BHRs is 
impossible. This gives an opportunity to consider these BHRs as one of the main 
candidates for cold dark matter in our Universe. 
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Figure 1. Illustration of the moment of last transition. Prelast state is characterized by the 
regular horizon with the usual quasi-Schwarzschild configuration. The last state is singular 
configuration making the transition from prelast to last state forbidden. 
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Figure 2. Shape of BH mass lost rate versus BH mass in Gauss-Bonnet case when the energy 
conservation is taken into account. Right part of the graph represents the usual Hawking 
evaporation law when —dM/dt ~ 1/M 2 . Left part shows the picture at last stages when 
evaporation decelerates and then stops, distinguishing the minimal possible mass ( "ground state" ) . 




Figure 3. Metric function a and e as a function of the mass M in Planck units for a fixed 
minimal mass Mum = WMpi. Stars are numerically computed values and the line is the fit used 
to derive the spectrum. 
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Figure 4. Im(S) (dots) and the fit (5 • 10~ 4 ) * (M — M mm ) 3 (dashed line) versus BH mass M 
during the last stages of BH evaporation in the Gauss-Bonnet case with M min = 10.6Mp;. It is 
necessary to note that during last stages of evaporation the emitted energy ut < M — M m i n <C 1 . 
For fixed values of u = u>* in the vicinity of M m i n (0(M m i n ) = 0.01) the mass M € (M m i n + 
ui*,M min +uj* +0(M mm )). So, for different values of u* = w* + 0(M min ),uj* 1 = 0.1, i e N) 

M belongs to different (without intersection) intervals. Finally, Im(S) is represented as connection 
of such intervals with the most probable values of w* G (0.1, 1.1). 



